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ABSTRACT 

While observations of large-scale structure and the cosmic microwave background (CMB) provide 
strong constraints on the amplitude of the primordial power spectrum (PPS) on scales larger than 
10 Mpc, the amplitude of the power spectrum on sub-galactic length scales is much more poorly 
constrained. We study early structure formation in a cosmological model with a blue-tilted PPS. We 
assume that the standard scale-invariant PPS is modified at small length scales as P{k) ~ with 
rus > 1. We run a series of cosmological hydrodynamic simulations to examine the dependence of 
the formation epoch and the characteristic mass of primordial stars on the tilt of the PPS. In models 
with rris > 1, star-forming gas clouds are formed at z > 100, when formation of hydrogen molecules 
is inefficient because the intense CMB radiation destroys chemical intermediates. Without efficient 
coolant, the gas clouds gravitationally contract while keeping a high temperature. The protostars 
formed in such “hot” clouds grow very rapidly by accretion to become extremely massive stars that may 
leave massive black holes with a few hundred solar-masses at z > 100. The shape of the PPS critically 
affects the properties and the formation epoch of the first generation of stars. Future experiments of 
the CMB polarization and the spectrum distortion may provide important information on the nature 
of the first stars and their formation epoch, and hence on the shape of the small-scale power spectrum. 
Subject headings: cosmology: theory - dark ages, reionization, primordial stars - stars: Population 
III - stars: formation - method: numerical 


1. INTRODUCTION 

The formation of the first stars marks the end of the 
cosmic Dark Ages, when their energetic photons initi¬ 
ates reionization of the inter-galactic medium. The ex¬ 
act epoch when cosmic reionization began is still uncer¬ 
tain, but evidence is mounting from an array of obser¬ 
vations that the first luminous objects su ch as galaxies 
and quasars were formed very early on (iWatso n et al.l 
201fiyWiet_^ 120151 : iPlanck Collaboration XIIJ 120151) . 
Wu et al.l (12015 1 recently reported the discovery of a su- 
permassive black hole in an extremely bright quasar. 
While their exact estimate of a black hole mass of twelve 
billion solar masses rests on the controversial use of broad 
emission lines (see e.g.. iCackett et all 1201511 . this ultra- 
luminous quasar must host a supermassive black hole. 
There are now over 40 quasars known with z > 6 in¬ 
cluding ULA S J1120-I-0641 with a 2 x 10^ black hole 
at z = 7.085 (jMortlock et al.l[2011l : iDe Rosa et al.ll2014[l 
and a recently discove red z = 6.889 quasar with a black 
hole of 2.1 X 10® Mq (|De Rosa et al.11201411 . Assembling 
these giant quasar hosts by this early epoch poses a se¬ 
rious challenge to the theory of structure formation and 
the growth of black holes. 

There has been significant progress in the past decade 
in theoretical studies of the format i on of the first gen- 
eration of stars (see iBromml 1201,81 iGloved 1201,81 : iGreifl 
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120151 for recent reviews). Owing largely to the fact 
that the initial conditions are cosmologically well deter¬ 
mined, one can perform ab initio simulations of early 
structure formation (e.g., iBromm et al.l[200^ lAbel et al.l 

l2002l : rYoshida et al.ll2003f) . An important element of such 

simulations is the primordial power spectrum (PPS) from 
which a realistic initial density field can be generated in 
a fully cosmological context. 

The PPS is observationally determined to have 
a power-law with a spectral index rig — 0.96 
(|Planck Collaboration XIll 1201511 . Unfortunately, such 
tight constraints are derived only at large length scales, 
up to the wavenumber of fc ~ 0.2 Mpc“^ (iHlozek et al.l 
1201211 . Most of theoretical studies thus assume a 
scale-invariant PPS and adopt significant extrapola¬ 
tion of the observationally determined large-scale power 
spectrum. The observational evidence on sub-galactic 
scales is confusing: milli-lensing measurements sug¬ 
gest substructure consistent with the extrapolation of 
the A-Cold Dark Matter (A-CDM) power spectrum 
to smaller scales (iHezaveh et al.l 120131 iBussmann et al.l 
[MlIxF et al. I [mi" yet studies of local dwarf galax¬ 
ies sug gest significant prob lems with A-CDM on small 
scales (|Strigari et al.l [200711 . It is not only interesting 
to explore how early structure formation is affected by 
enhanced density fluctuations at small scales, but also 
timely because the existence of supermassive black holes 
at very high redshifts suggests rapid growth of small scale 
structure in the early universe. 

A variety of possibilities are proposed from the 
physics of the early universe that posit scale-dependent 
PPS. Popular single-field inflation models predict a 
nearly scale-free and slightly red power spectrum 
for adiabatic scalar perturbations. There are also 
variant models that predict scale-dependent power 
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Fig. 1. — We plot the matter power spectra for the considered 
cases in the present paper with kp = 100 h Mpc“^. The black line 
is for the scale-independent PPS (Eq. 0 whereas the other three 
lines show the scale-dependent models given by Eq. [2] with rris = 
1.5 (blue), 2.0 (green), and 2.5 (red), respectively. 


spectrum with kinks, bumps, a nd other typ es _ of 

features fe.g.. iStarobins^ Il992t lAdams et alJ 119971: 

Kawaguchi et al.l I2008I iBiswas et al. I l201Ct iBarnabvl 

2010li . Of particular interest to us in the present pa¬ 
per are the models that generate density perturba¬ 
tions with a blue-tilted power spe ctrum or with en¬ 
hancement at small length scales (iCovi fc Lyth I 119991 : 
iMartin fc Brandenbcrgedl200lHGong fc Sasakill2011h ~ 

It is naively expected that structure forms early in 
such models, yielding abundant small mass dark halos. 
Because primordial star formation involves a number of 
physical processes such as chemistry and radiative trans¬ 
fer as well as gravitational assembly of dark halos, the 
initial density perturbations or the PPS can affect the 
formation of primordial stars in a complicated manner. 
Clearly, it is important to study the properties of the 
first structure using cosmological simulations with the 
relevant physics included. 

In the present paper, we perform a series of cosmo¬ 
logical simulations for different PPS models to study 
early structure formation in detail. In particular, we 
examine how the slope of the PPS at the small length 
scales affects the formation epoch and the mass of the 
first generation of stars. We begin by describing the 
calculation methods in Section [2l Section |3] shows the 
simulation results of primordial star formation for differ¬ 
ent PPS models. We discuss the dependence of the pri¬ 
mordial star formation on the PPS model in Section 01 
Throughout this paper, we adopt the standard A-CDM 
cosmology with the total matter density Ilm = 0.3086, 
the baryon density fib = 0.04825, the dark energy den¬ 
sity Ha = 0.6914 in units of the critical density, the 
Hubble cons tant h = 0.6777, and the primo rdial index 
rig = 0.9611 (|Planck Collaboration XVIII2014D . 

2. NUMERICAL SIMULATIONS 

We perform a series of cosmological simulations that 
start from the initial conditions generated for differ¬ 
ent PPS models Pprim(fc)- In the next sections, we 



Fig. 2. — We compare the halo mass functions a.t z = 100 for 
three PPS models. We set the normalization parameter as = 2.5 
here. The lines show the Press-Schechter mass functions, while 
the symbols show the results of our cosmological simulations with 
a volume of 60 h~^ comoving kilo-parsec on a side. We run the 
Friend-Of Friend (FOF) halo finder with linking length b = 0.2 to 
obtain the halo mass functions. 


describe the numerical methods with emphasis on the 
improvement over our earlier numerical simulations 
(|Hirano et al.l[2?)T^ l2015[ ). 

2.1. Initial Conditions 

The standard scale-invariant PPS is given by 

T’prim(fc) c>cfc"= , (1) 

whereas the model power spectrum with enhancement at 
small scales is given by 

Pprim(fc)oc V" (fc < fcp) , 

ocfc”^-’”^ •fc’”^ (fc > fcp) . (2) 

We adopt two parameters that characterize the spec¬ 
trum; fcp is the pivot wavenumber above which the spec¬ 
tral index is set to be Wg. We consider fcp = 100, 300, 
and 500 h Mpc“^ and rus = 1.5, 2, and 2.5. 

Figure [T] shows the matter power spectra for models 
with fcp = 100 h Mpc“^. The fiducial pure power-law 
PPS case is also show n there for comparison. We use 
the transfer function of lELsenstein fc Hul (jl999f ) for this 
plot and also for generating the initial conditions us¬ 
ing t he publicly available code MUSIC (jHahn fc Abell 
1201 If) . The density fluctuation amplitude is normalized 
by (Ts- Because our simulation volume is very small, we 
set CTs larger than the current stan dard value of about 0.8 
(jPlanck Collaboration XIIIIl2015h . so that we can simu¬ 
late the formation of high-cr peak halos. In practice, we 
set (Tg = 1.5, 2, and 2.5. In total, we generate 18 cosmo¬ 
logical initial conditions with the parameter sets given in 
Table [H 

The effect of the enhanced small-scale power is clearly 
seen in the mass function (Fig. [2]). The halo mass func¬ 
tions differ by orders of magnitude at z = 100. Fig¬ 
ure 0] also shows that the Press-Schechter mass function 
agrees fairly wel l with the simulation results (see also 
iReed et alH2005D . 
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TABLE 1 

Properties of First Star Formation Depending on the PPS Models 


ID 

0-8 

Alp 

[h Mpc“^) 

mg 

•^form 

Mvir 

(10® Mq) 

-^Jeans 

(Mq) 

-Mjeans 

(10-3 Mq yr-1) 

Mf^caX 

(Mq) 

(Mq) 

too 

1.5 

- 

- 

24.0 

2.17 

45 

0.63 

45 

88 

111 


100 

1.5 

27.1 

1.82 

162 

2.97 


261 

112 



2.0 

36.1 

0.88 

162 

1.16 


135 

113 



2.5 

87.9 

0.31 

69 

0.77 

96 

101 

131 


300 

1.5 

24.0 

2.11 

23 

0.24 


45 

132 



2.0 

25.2 

2.19 

124 

0.74 


98 

133 



2.5 

33.9 

1.97 

171 

3.31 


281 

151 


500 

1.5 

23.5 

2.10 

33 

0.35 


58 

152 



2.0 

24.0 

2.19 

61 

1.02 


123 

153 



2.5 

27.1 

2.83 

90 

1.61 


170 

200 

2.0 

— 

_ 

34.9 

1.30 

401 

1.54 

114 

165 

211 


100 

1.5 

40.1 

0.83 

167 

1.29 


145 

212 



2.0 

53.2 

0.22 

57 

0.25 


46 

213 



2.5 

138.2 

0.24 

142 

1.71 

145 

177 

300 

2.5 

_ 

_ 

46.2 

0.83 

243 

2.34 

118 

220 

311 


100 

1.5 

52.6 

0.65 

164 

1.39 


153 

312 



2.0 

72.1 

0.22 

60 

0.43 


67 

313 



2.5 

186.3 

0.37 

411 

2.72 

362 

245 


Note. — Column 1: Cloud index, Column 2-4: PPS model parameters, Column 5-8: Properties of star-forming clouds. Column 9: 
Stellar mass calculated by RHD simulations, and Column 10: Stellar mass estimated by a fitting formula The cloud ID indicates the 
combination of the three parameters {cjs, kp, mg}. 


2.2. Cosmological Simulations 

We perform the cosmological simulations using the par¬ 
allel N-body / S moothed Particl e Hydrodynamics (SPH) 
code Gadget-2 (|Springell 1200511 suitably modified for 
the primordial star formation simula tions (jYoshida et al.l 
120061 l2007t iHirano et al.ll20l41201511 . In order to achieve 
sufficient numerical resolution, we use hierarchical zoom 
initial conditions. The parent simulation volume is 
Tbox = 100 h~^ kpc on a side. The mass of the refined 
dark matter particles in the zoomed region is 0.80 Mq, 
and that of the baryonic component is 0.15 Mq. We fol¬ 
low structure formation from Zini = 499 until the central 
hydrogen number density reaches nn.cen — 10^ cm“^ or 
10^^ cm ~^ in a gas cloud that first collapses (see Sec¬ 
tion [231). The init ial ionization fract i on is c omputed by 
using RECFAST (iSeaeer et al.l 119991 l2000t iWong et al.l 
(200^ as Xe = 6.88 X 10 ^ at the initial redshift. Dur¬ 
ing the cloud collapse, we keep the refinement criterion 
that the local Jeans length is always resolved by 15 times 
the local smoothing length by progressively increasing 
the spatial resolution by the p article-splitting technique 
(|Kitsionas fc Whitworth|l2002fl . 

Let us briefly describe the overall trend found from 
the simulation result. In the standard model, the first 
stars are formed at z 20 — 50 le.g.. iTegmark et al.l 
Il997t iReed et 31120051 : lYoshida et al.l I2003D . liT one of 
our simulations with as = 1-5, the first star is formed at 
-Zform = 24 within the simulation volume of 100 h~^ kpc 
on a side. We label this sample as ID 100. Here, the first 
digit “1” indicates the value of as (see Tabled]). We also 
find that, by changing as to 2.0 and 2.5 (ID = 200 and 
300), the first star in the simulation volume is formed 
earlier at Zform = 35 and 46, respectively. 

2.2.1. H 2 Formation in the Early Universe 

We are interested in the thermal evolution of a pri¬ 
mordial gas cloud formed at z > 100. Because CMB 


photons at such high redshifts are capable of dissociat¬ 
ing molecules and ions, several chemical reactions need to 
be included in order to follow the thermal and chemical 
evolution accurately. 

At z < 100, the main H 2 formation path is the so-called 
H“ channel; 

R + (3) 

H-+H^H2+e-. (4) 

Because H~ ions are destroyed by CMB photons at z > 
100 via the photo-detachment reaction 

R- +hiy^R + e- , (5) 

the H“ channel is not effective at such early epochs. An 
alternative path is the Hj channel: 

R+ + R^R+ + hv , ( 6 ) 

H++H^H2+H+. (7) 

An opposing dissociative process by energetic photons 
operates as 

R+ + hiy^R+ + R . ( 8 ) 

At z > 120, most of the hydrogen molecules are formed 
by the HJ channel. The mean H 2 fraction in the Universe 
reaches ~ 10“® at z ~ 60. We update the chemical 
ne twork by introducin g the above reactions (see a review 
bviGalli fc Pallall2013ll . The reaction rates are based on 
iCoDDola et al.l (|2013[) . We describe more details of the 
chemistry implementation in Appendix [^ 

2.3. Final Stellar Masses 

It is important to study the characteristic mass of 
the first stars in our variant cosmological models with 
blue-tilted PPS. We calculate the mass of a primordial 
star by using two different methods. One is direct cal¬ 
culation of gas accretion onto the central protostar us- 
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Pdm [cm-T Pgas [cm ’] [K] 


Fig. 3.— Projected densities for dark matter (Panel a) and baryonic components (b). We also plot the projected gas temperature (c) 
at the formation sites of the primordial stars ID 300 (top; Zf^rm = 46) and ID 313 (bottom; 2:form = 186). The plotted region has a proper 
100 pc on a side. We scale the density in units of particle number density per cubic-centimeter. 


ing the two-dimensional axisymmetric radiation hydro- 
dynami c (RHD~) simulation coupled with the stellar evo¬ 
lution (jHosokawa et alJ 120111 l2012t iHirano et alJ 12014 
1201511 . We perform the calculations for 6 cases. We 
label them as ID = 100, 113, 200, 213, 300, and 313. 
Here, the second digit of the object ID indicates the pivot 
wavenumber, and the last digit indicates the value of ms 
(see Tabled]). For each gas cloud, we generate the initial 
conditions for the RHD simulation by suitably mapping 
the final snapshot of our cosmological simulation at the 
time when the density of the collapsing cloud core reaches 
?T-H.cen = 10^^ cm“^. The RHD simulations are run until 
the mass accretion rate onto the central star falls below 
10“^ Mq yr“^. At this point the stellar growth is al¬ 
most completely halted and the resultant mass can be 
regarded as the final stellar mass. We denote the stel¬ 
lar masses calculated directly by the RHD simulation as 

We also estimate the final stellar masses by using a 
fitting function that is derived from a number of RHD 
simulations of accreting protostars in our earlier stud¬ 
ies. The fitting function essentially_ correlates the stellar 
mass M*_est to the gas infall rate Mjeans evaluated for a 
gravitational l y Jean s unstable gas cloud. It is given by 
IHirano et al.l ()20I5[1 as 


M*,est = 250 


I Jeans 

I 2.8 X 10-3 Mq yr- 


By adopting this formula, we can skip the numerically 
costly RHD simulation in the accretion phase. Note that 
the difference between the calculated and estimated val¬ 
ues is about fifty percent, and mostly within a factor of 


two (see £g. 14b in IHirano et al.l [20141 . Although this 
might still appear insufficient to determine the stellar 
mass for all the cases in Table dl the mass variations 
caused by model parameters (kp,ms, etc) are typically 
larger than the uncertainty in the estimated mass (see 
Table dj • Because we aim at investigating how cosmo¬ 
logical parameters such as the slope of the power spec¬ 
trum affect the formation of early halos and gas clouds, it 
would be clearer to interpret the resulting stellar masses 
in terms of gas clouds’ properties such as the Jeans mass 
and accretion rate. Therefore, we use the estimated stel¬ 
lar masses (Eq. [9]) for all the 18 cases to discuss the 
qualitative dependence on the adopted PPS models. 

Our stellar mass estimates are based on the results 
of radiation-hydrodynamics simulations which follow 
only the evolution of the central star. Recent three- 
dimensional simulations show that the circumstellar disk 
can gravitationally f ragment to produce multiple clumps 
(sink particles: e.g..Oark et al.ll2011l : IGreif et HI 1201 IL 
120121 : IStacv fc Brommll2013ll . Such fragments, if some of 
them survive over a long period of disk accretion, would 
form multiple stars with correspondingly small masses. 
High-resolution simulations show also that fragments ef- 
hciently merge onto the central star in roughly an or¬ 
bital timescale (e.g., IGreif et HI 120121 : IVorobvov et al.l 
I2013f) . In the latter case, the final stellar mass would 
be similar with the calculated value by our radiation- 
hydrodynamics simulation. Gonservatively, we interpret 
the derived M* as the total mass of stars formed, rather 
than the mass of a single star. 

3. RESULTS 
3.1. Major Effects 
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Redshift Redshift 


Fig. 4.— Time evolution of gas density (Panel a) and H 2 molecule fraction (Panel b) of the densest gas elements for cases with ag = 2.5 
until nH,cen = 10® cm~®. The lines with different colors represent the four PPS models in Table [T] 


We find major effects of the small-scale power spec¬ 
trum on the properties of host halos and on the star for¬ 
mation within them. Figure [3] shows the local conditions 
around the formation site of the primordial star. The top 
panels are for the case with normal scale-invariant PPS 
(ID 300). The primordial star-forming cloud is formed 
at 2form = 46 inside a dark matter mini-halo that has 
a virial radius of i?vir — 19 pc. The bottom panels are 
for the case with an extremely enhanced PPS (ID 313). 
The gas cloud is formed earlier at Zform = 186 inside a 
more compact mini-halo with i?vir — 3.6 pc. Figure[3Ja) 
shows the density distribution of dark matter compo¬ 
nent within a proper 100 pc on a side. Note that the 
mean density is higher in the bottom case than in the 
top because of the cosmic expansion, as p oc (1-1- z)^. 
The number of small mass dark matter clumps is signifi¬ 
cantly larger in the bottom panel, reflecting the enhanced 
initial density perturbations at the corresponding length 
scales (e.g.. Figs. [1] and [2|). Interestingly, we find a sec¬ 
ondary star-forming cloud formed in the vicinity of the 
first one within several 10 pc (Fig. [SJd). Such a close 
pair of primordial stars likely affect each other through 
radiative, dynamical, and chemical feedback effects. In 
the top panels, we also find a neighboring mini-halo but 
it has not grown sufficiently to host a dense star-forming 
cloud (see also the temperature distribution in Panel c). 

The host mini-halos forming at very high redshifts are 
generally more compact and have larger dark matter 
densities, whereas the characteristic density of the gas 
cloud is determined largely by molecular hydrogen cool¬ 
ing. Consequently, the baryonic component becomes dy¬ 
namically dominant only after the gas density exceeds 
~ 10^ cm“^ in very early halos. 

3.2. Formation Epoch and H 2 Formation 

Figure SJa) shows the density evolution of the gas 
cloud’s core as a function of redshift for the cases with 
(Tg = 2.5. With more enhanced power (larger Wg), 
the clouds collapse earlier. In the extreme case with 
ms = 2.5 (red line), the cloud collapses at 2 : 186. Then 

H 2 formation mainly proceeds via the HJ channel, but 


the photo-dissociation by CMB prevents H 2 formation. 
Hence the molecular fraction remains low for ID 313 until 
three-body H 2 formation operates at very high densities 
(Fig. Hb). The formation of H 2 molecules critically af¬ 
fects the thermal evolution, prestellar collapse, and the 
subsequent accretion process onto the central protostar. 

3.3. Thermal and Dynamical Properties of Clouds 

Figure [S] shows the evolution and the properties of col¬ 
lapsing clouds for the six cases for which direct RHD 
simulations are performed. We use black lines for mod¬ 
els with the standard power-law spectrum (mg = rig) and 
we use red lines for models with enhanced small scale 
power (mg > 1). The thermal evolution in the standard 
model is well understood and can be described as follows 
(see Panel a). The temperature first increases due to 
gravitational contraction but then decreases owing to H 2 
cooling, when the molecular fraction reaches the critical 
value that is needed to cool the ga s within a Hubble tirn e 
/h 2 — 5 X 10“"^ (see Fig. [SJ), e.g.. iTeemark et aI1ll997tl . 
Note that the minimum temperature is floored by Tqmb 
at high redshifts. The gas temperature starts increas¬ 
ing again when the compressional heating overcomes the 
radiative cooling, leading to the formation of a protostar. 

In our models with enhanced small scale power (red 
lines), the cloud evolves on a higher temperature track 
because of the lower H 2 fraction (Panel b). With larger 
mg, the density at which the H 2 fraction reaches the 
critical value for cooling is larger. Interestingly, for ID 
313 (red solid line), the H 2 fraction does not reach the 
critical value until three-body H 2 formation occurs at 
■'T-H.cen > 10® cm“®. In this case, the temperature at the 
onset of run-away collapse is 700 K, much higher than 
that found in the usual primordial star formation. One 
can naively expect that protostars formed in such a ’hot’ 
cloud will grow very rapidly. 

The stellar mass can be estimated from Eq. ([9|) . Effec¬ 
tively, the mass accretion history determines the proto- 
stellar evolution and the strength of radiative feedback 
that halts the accretion of surrounding medium. Fig¬ 
ure [SJc) shows the gas infall rate as a function of the 
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nH,cen Mgn^(r) [Mg] 


Fig. 5. — The evolution of temperature (Panel a) and H2 fraction (Panel b) at the center of the cloud. The horizontal lines in Panel (a) 
indicate the CMB temperature, Tcmb(-z) = 2.73 (1 + z), at the respective formation epoch ^form listed in Table[T] The dotted lines show 
the Jeans mass for 10 — 10® Mq in the density-temperature plane. Panels (c) and (d) display the profiles of the instantaneous gas infall 
rate and the degree of rotational support (/Kepler = I'rot/uKepler) 3-s a function of the enclosed mass. The profiles are calculated when the 
cloud core density reaches njj cen = 10^® cm“*. 


enclosed gas mass. On average, the gas infall rate sys¬ 
tematically increases at higher redshift (see the difference 
from long-dashed to solid lines). However, the gas infall 
rate does not monotonically increase with ms- This is 
because the host mini-halos of the gas clouds are formed 
at earlier epochs but contain less amount of gas. 

Figure EKd) shows the ratio of rotational speed to the 
Keplerian speed. Interestingly the gas clouds formed at 
earlier epochs spin significantly more slowly in the mod¬ 
els with enhanced power. This may reflect the effective 
redistribution of the (orbital) angular momenta brought 
in by numerous minor mergers. Stars formed in a gas 
cloud with less spin tend to be more massive because 
of relatively rapid s pherical accretion onto the protostar 
(|Hirano et al.ll20f^ . 

3.4. Dependence on Model Parameters 

It is important to study the effect of the pivot scale 
kp of the PPS on gas cloud formation. There is an im¬ 
portant physical scale that characterizes primordial star 
formation. Let us denote the characteristic wavenumber 
^haio for a mini-halo that has a physical size i?haio so 
that fchaio ~ 27r/i?haio- When kp is smaller than fchaio, 
the enhanced power generates more structure of mini¬ 
halo sizes. 

In the fiducial case of ID = 100, the host halo has 
a virial mass of ~ 2 x 10® Mq. The corresponding 
wavenumber is fchaio = 27r/i?haio ~ 500 h Mpc“^. When 
kp is comparable to fchaio, as in the cases with kp = 300 
and 500 h Mpc“^, increasing mg does not significantly 
affect 2:form and Mvir (see Table [T|). This is in direct con¬ 
trast to the trend found when kp < fchaio! increasing rrig 


actually resulted in larger Zform and smaller Mvir. This 
can be understood as follows. The formation epoch is 
largely determined by the amplitude of large-scale den¬ 
sity perturbations whereas the enhanced power promotes 
formation of small mass halos. Therefore, when kp is 
small, density perturbations at all the relevant length 
scales are enhanced. The mini-halos are then formed 
early, and are physically compact. 

It is worth considering possible halo-to-halo variation. 
Our underst anding so far is based on a limited number 
of samples. iHirano et al.l (|20I4 I20I5I1 show substantial 
differences in the properties of star-forming clouds even 
with similar masses and similar formation epochs that 
are in turn determined largely by the cosmological pa¬ 
rameters. In principle, using a much large number of 
halos could show the general trend and variations more 
clearly, but we have focused on qualitative differences 
caused by the shape of PPS. 

4. DISCUSSION 

We have studied how the shape of the primordial power 
spectrum at small scales affects early structure forma¬ 
tion. For a significantly enhanced PPS, star-forming gas 
clouds are formed at extremely early epochs of z > 100 
(ID 213 and 313), when the usual gas-phase reaction 
of molecular hydrogen formation does not operate. At 
100 < z < 300, hydrogen molecules are formed via the 
channel and thus the gas can still cool and condense 
to form primordial stars. The formed stars have masses 
as large as 300 Mq. 

At even earlier epochs of z > 300, there is no effi¬ 
cient formation channel of H 2 and the possible path to 
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star formation is through cooling by atomic hydrogen. If 
the highly enhanced density perturbations, as explored 
in the present paper, grow quickly to form primordial 
gas clouds in the very early universe, an extremely mas¬ 
sive star might be formed in a similar manner to the so- 
called direct collaps e black hole formation model (e.g., 
Ilnavoshi et al.l[20l3l . Such very massive stars gravita¬ 
tionally collapse to form remnant massive black hole at 
the end of their lives, leaving promising “seeds” for the 
formation of the recently dis covered supermassive black 

hole (SMBH;|meEaD[20il)- 


There is an additional generic physical process that can 
significantly affect the structure formation in the early 
universe; the relative motions between dark matter and 
baryons caused by acoustic oscillations at the cosmic re¬ 
combination era ( Tsellakhoyich^^^kat^IlQin) . The rel¬ 
ative streaming motions prevent the gas contraction into 
dark matter mini-halo and delay the star formation, and 
thus possibly chan ge the physic al condition of the first 
star formation fsee lFialkovll20l4 for recent review). The 
effect is important at the very high redshifts considered in 
this paper, since the root-mean-square streaming veloc¬ 
ity scales with redshift as oc (1 -I- z). We have run several 
additional simulations and confirmed that, under signifi¬ 
cant streaming motions, primordial gas clouds are formed 
later in large host halos than presented in the present 
paper. There is an interesting possibility that stream¬ 
ing motions delay the first star formation to the extent 
that gas collapse is suppressed in mini-halos and more 
massive halos eventually host very massive stars, which 
can become the seed of SMBH (iTanaka fc Lill20f4) . The 
blue-tilted primordial power-spectral can effectively in¬ 
crease the number of massive halos in the early universe. 
Further results from simulations with the streaming mo¬ 
tions will be presented elsewhere (Hirano et ah, in prepa¬ 
ration) . 

Observationally, an independent probe of the PPS at 
small length scales might be provided by the existence of 
small mass dark halos i n and around galax ies. Following 
the analytical model of iZhao et al.l (1200511 . we estimate 
that an extremely large number of halos with masses of 
1 Mq are formed for uis = 1.5 and kp = 100 Mpc“^. 
The typical formation epoch of such solar-mass halos is 
z ^ 300, and thus the halos may remain as very compact 
objects. However, even the dense and compact clumps 
can be destroyed by tidal stripping/disruption and en¬ 


counter with stars in the Galaxy. It is thus difficult to 
estimate the abundance of the mini- and micro-halos sur- 
viving until today, either by large scale simula tion (e.g., 
iDiemand et al.l[20^ iGree n fc Go odwinl 120071 1 or by an¬ 
alytical calculation (IZhao et al.l 1200511 Nevertheless, if 
signatures of small mass dark halos are detected through, 
for instance, dark matter annihilation or direct detection 
experiments, one can obtain invaluable information on 
their origin, i.e., the small-scale power spectrum. 

Currently, there are little direct probes of the small- 
scale primordial density fluctuatio ns, and only model 
dependent constraints are available (iCh uba et al.ll2012l : 
iBerezinskv et al.ll20l4 iNataraian et al.l 201511 . We look 
forward to future CMB expe riments such as PRIS1V0 
and PIXIE (|Kogut et al.ll201lll that can infer the forma¬ 
tion epoch and the spectral energy distribution (hence 
roughly the typical mass) of the first stars through mea¬ 
surement of CMB spectral distortions and the detec- 
tion of low level o f early reionization at Z > 10 (e.g., 
iRicotti et al.l 1200511. Finally . James Webb Space Tele- 
scone (iWindhorst et al.l200^ and WFIRST observations 
(|Spergel et al.ll2015ll will probe the existence of quasars 
at even higher redshifts. Understanding the nature of 
the first stars may thus provide a route to probe obser¬ 
vationally the small-scale primordial power spectrum. 
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APPENDIX 

A. REACTION RATE CALCULATION OF H+ PHOTO-DISSOCIATION 

The react ion rate calculation of photo-d issociation is h i ghly u ncertain. In our previous calculations (iHirano et alJ 
we adopt the fitting function of iGalli fc Fallal (ii998f) assu ming the LTE-level population . The accurate 
react ion rate can be calculated by considering the HJ-level population (iHirata fc Padnianabhanl[2006t iCoppola et alJ 
1201 It) . The time evolution of H 2 abundance at z > 60. which is the releva nt epoch in our study, differs significantly 
depending on the adopted reaction rate (see fig. 3a in iGalli fc Pallal[2013l ). We take a straightforward approach by 
comparing simulations with different reaction rates. Fortunately, we find that the uncertainty in the reaction rate 
gives essentially no influence on the primordial star formation. Stars are formed in dense gas clouds in which a critical 
amount of hydrogen molecules are formed to overcome the compressional heating by their radiative cooling. When a 
gas cloud becomes gravitationally unstable after the so-called loitering phase, the H 2 fraction within the cloud reaches 
a threshold value that is independent on the treatment of photo-dissociation. We thus adopt the LTE rate to 
produce the main results in the present paper. In most cases, we find a constant “floor” of the H 2 fraction at 2 x 10“^ 
owing to the H 2 formation via the channel (Fig. IHd). 



















